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ABSTRACT 

We constructed a sample of 185 Flat Spectrum Radio Quasars (FSRQs) by cross- 
correlating the Shen et al.'s SDSS DR3 X-ray quasar sample with FIRST and GB6 
radio catalogues. From the spectrum energy distribution (SED) constructed using 
multi-band (radio, UV, optical, Infrared and X-ray) data, we derived the synchrotron 
peak frequency and peak luminosity. The black hole mass A/bh and the broad line 
region (BLR) luminosity (then the bolometric luminosity iboi) were obtained by mea- 
suring the line-width and strength of broad emission lines from SDSS spectra. We 
define a subsample of 118 FSRQs, of which the nonthermal jet emission is thought to 
be dominated over the thermal emission from accretion disk and host galaxy. For this 
subsample, we found 25 FSRQs having synchrotron peak frequency v^cak > 10^^ Hz, 
which is higher than the typical value for FSRQs. These sources with high I'poak could 
be the targets for the Fermi Gamma-ray telescope. Only a weak anti-correlation is 
found between the synchrotron peak frequency and peak luminosity, while no strong 
correlation is present either between the synchrotron peak frequency and black hole 
mass, or between the synchrotron peak frequency and the Eddington ratio iboi/^Edd- 
When combining the FSRQs subsample with the Wu et al.'s sample of 170 BL Lac 
objects, the strong anti-correlation between the synchrotron peak frequency and lumi- 
nosity apparently presents covering about seven order of magnitude in i^poak- However, 
the anti-correlation differs with the blazar sequence in the large scatter. At similar peak 
frequency, the peak luminosity of FSRQs with v-pcak > 10^^ Hz is systematically higher 
than that of BL Lac objects, with some FSRQs out of the range covered by BL Lac 
objects. Although high Vpcak are found in some FSRQs, they do not reach the extreme 
value of BL Lacs. For the subsample of 118 FSRQs, we found significant correlations 
between the peak luminosity and black hole mass, the Eddington ratio, and the BLR 
luminosity, indicating that the jet physics may be tightly related with the accretion 
process. 

Key words: galaxies: active — galaxies: quasars — quasars: emission lines — quasars: 
jet 



1 INTRODUCTION 

Blazars, including BL Lac objects and flat-spectrum radio quasars (FSRQs), are the most extreme class of active galactic 
nuclei (AGNs), characterized by strong and rapid variability, high polarization, and apparent superluminal motion. These 
extreme properties are generally interpreted as a consequence of non-thermal emission from a relativistic jet oriented close to 
the line of sight. As such, they represent a fortuitous natural laboratory with which to study the physical properties of jets, 
and, ultimately, the mechanisms of energy extraction from the central supermassive black holes. 



* E-mail:zychen@shao. ac.cn 



2 Z. Y. Chen, M. F. Gu and X. Cao 



The most prominent characteristic of the overall spectral energy distribution (SED) of blazars is the double-peak structure 
with two broad spectral components. The first, lower frequency component is generally interpreted as being due to synchrotron 
emission, and the second, higher frequency one as being due to inverse Compton emission. BL Lac objects (BL Lacs) usually 
have no or only very weak emission lines, but have a strong highly variable and polarized non-thermal continuum emission 
ranging from radio to 7-ray band, and their jets have synchrotron peak frequencies ranging from IR/optical to UV/soft-X-ray 
energies. Compared to BL Lacs, FSRQs have strong narrow and broad emission lines, however generally have low synchrotron 
peak frequency. According to the synchrotron peak frequency, BL Lac objects can be divided into three subclasses, i.e. low 
frequency peaked BL Lac o bjects (LBL), intermediate objects (IBL) and high frequency peaked BL Lac objects (HBL) 
JPadova ni fc Gi ommi 19951 ). I n general, radio-selecte d BL Lacs tend to be LBLs, and XBLs are HBLs ( Urrv fc Padovani 



199 



I. 



Fossati et al 



1 19981 ) and Ghisellini et al. 1 19981 ) have proposed the well-known 'blazar sequence' , which plot 



various 

powers vs the synchrotron peak frequency I'peak including FSRQs and BL Lacs. They draw a conclusion that the peak frequency 
seem to be anti-correlated with the source power with the most powerful sources having the relatively small !^peak and the 
least powerful ones having the highest I'peak- Ghisellini et al. ( 19981 ) gave a theoretical interpretation to these anti-correlations, 
namely, the more powerful sources suffered a larger probability of losing energy, the more cooling the sources subjected, thus 
translates into a lower value of t'peak. 

However, recently the blazar sequence has been largely in debates. By constructing a large sample of about 500 blazars from 



the D eep X-Ray Radio Blazar Survey (DXRBS) and the ROSAT All-Sky Survey-Green Bank Survey fRGB). IPadovani et al 



( 20031 ) found that the "X-ray-strong" radio quasars, with similar SED to that of HBLs, have much higher synchrotron peak 
frequencies than those of classical FSR Qs. Their DXRBS s ample does not show the expected blazar sequence. Exceptions 
to blazar sequences were also found by lAnton fc Browne I (|2005l ) that most of the low radio luminosity sources have syn- 
chrotron peaks at low frequencies, instead of the expected high freq uencies in their blazar sample. They claimed that at 



least part of the systematic trend seen by Fossati et al. ( 19981 ) and Ghisellini et al. ( 19981 ) results from selection effects 



Nieppola. Tornikoski fc Valtaoia ( 20061 ) have studied a sample of over 300 BL Lacs objects, of which 22 objects have high 



i-peak > 10^® Hz. There are negative correlations between i/peak and the luminosity at 5 GHz, 37 GHz, and 5500A, however, 
the correlation turns to slightly positive in X-ray band. Moreover, they claimed that there is no significant correlation between 
source luminosity at synchrotron peak frequency and t^peak, and several low energy peaked BL Lacs with low radio luminosity 
were also found. By using the results of very recent survevs. lPadovani I (|2007l ) claimed that there is no anti-correlation between 
the radio power and synchrotron peak frequency in blazars once selection effects are properly taken into account, and some 
blazars were found to have low power as well as low i/poak, or high power an d high t'peak as well. Furthermore, FSRQs with 
synchrotron peak frequency in the UV/X-ray band have been claimed (e.g. Giommi et al.l 2007 : Padovani et al. 20021) . In 



summary, it seems that the blazar sequence in its simplest form cannot be valid ([Padovani 



I2OO 



In this paper, we investigate the dependence of synchrotron peak frequency on the peak luminosity for a sample of FSRQs 
selected from SDSS DR3 quasar catalogue. The advantage of our sample is that the SDSS spectra enable us to measure the 
various broad emission lines, and then to estimate the black hole mass and the broad line region luminosity. Therefore, the 
relationship between the synchrotron peak frequency and black hole mass and Eddington ratio can be explored, which might 
have potential importance for us to study the jet formation and physics and the jet-disk relation. We present the sample 
selection in § 2. The reduction of SDSS spectra and the estimation of black hole mass and broad line region luminosity are 
described in § 3. The derivation of synchrotron peak frequency and luminosity are given in § 4, in which the thermal emission 
from accretion disc and host galaxy are also calculated. The various correlation analysis are shown in § 5, of which we focus on 
the relationship between the synchrotron peak frequency and luminosity. § 6 is dedicated to discussions. Finally, the summary 
is given in § 7. The cosmological parameters Ho = 70km s~^ Mpc~^, f2m=0.3, S7a = 0.7 are used throughout the paper, and 
the spectral index a is defined as oc v^" with being the fiux density at frequency v. 



2 THE SAMPLE SELECTION 



2.1 The SDSS DR3 X-ray quasar sample 

We started from the SDSS DR3 X-ray quasar sample of Shen et al. ( 20061 ) . which is the result of individual X-ray detections of 
SDSS DR3 quasar catalogue (Schneider et al. 2005) in the images of ROSAT All Sky Survey (RASS). The SDSS DR3 quasar 
catalogue consists of 46,420 objects with luminosities brighter than Mi = —22, with at least one emission line with full width 
at half-maximum (FWHM) larger than 1000 km s~^ and with highly reliable redshifts. A few unambiguous broad absorption 
line quasars are also included. The sky coverage of the sample is about 4188 deg'^ and the redshifts range from 0.08 to 5.41. The 
five-band (u, g, r, i, z) magnitudes have typical errors of about 0.03 mag. The spectra cover the wavelength range from 3800 
to 9200 A with a resolution of abou t 1800 - 200 (see Schneider et al. 2005 for details). The soft X-ray prope rties of SDSS DR3 
quasars have been investigated by Shen et al. ( 2006[ ) through individual detection and stacking analysis. Shen et al. ( 20061 ) 
applied the upper-limit maximum likelihood method to detect the X-ray flux at the position of each SDSS DR3 quasar and 
accept the objects with detection liklihood L > 7 as individual detections. The number of these individual X-ray detected 
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quasars were 3366, which is about 25 percent higher than the RASS catalogue matches (see Ishen et alj [20061 . for details). 
The 1 keV X-ray luminosity in the source rest frame of this 3366 X-ray quasar sample is obtained by assuming a power-law 
distribution of X-ray photo ns N(E) oc E"*" wit h F ^ 2 and corrected for absorption u sing the fixed column density at the 
Galactic value according to Dickev fc Lockman ( 1990l l for each source (|Shen et al.ll2006l ). 



2.2 Cross-correlation with FIRST and GB6 radio catalogues 



In this paper, we define a quasar to be FSRQ according to the radio spectral index. Therefore, we cross-correlate the SDSS 
DR3 X-ray quasar sample wit h Faint Images of the Radio Sky at Twenty-Centimeters 1.4 GHz radio cata logue (FIRST, 
Becker .White fc Helfandlll995l ) and the Green Bank 6 cm radio survey at 4.85 GHz radio catalogue (GB6, ICregorv et al.l 
19961 ) ■ which are two of the largest radio surveys well matched with SDSS sky coverage. The FIRST survey used the VLA to 
observe the sky at 20 cm (1.4 GHz) with a beam size of 5. 4. FIRST was designed to cover the same region of the sky as the 
SDSS, and it observed 9000 deg^ at the north Galactic cap and a smaller ~ 2.5° wide strip along the celestial equator. It is 
95% complete to 2 mjy and 80% complete to the survey limit of 1 mjy. The survey contains over 800,000 unique sources, with 
an astrometric uncertainty of < 1 . Due to the deeper survey limit and higher resolution, we prefer to use FIRST, instead of 
NVSS, which was also carried out using the VLA at 1.4 GHz to survey the entire sky north of 5 = —40° and contains over 
1.8 million unique detections brighter than 2.5 mJy, however with lower spatial resolution 45 beam~^. The GB6 survey at 
4.85 GHz was executed with the 91 m Green Bank telescope in 1986 November and 1987 October. Data from both epochs 
were assembled into a survey covering the 0° < 5 < 75° sky down to a limiting fiux of 18 mJy, with 3.5 resolution. GB6 
contains over 75,000 sou rces, and has a positional uncertainty of about 10 at the bright end and about 50 for faint sources 
( KimbaU fc Ivezicll2008h . 

The sample of 3366 quasars was firstl y cross-correlate d between the SDSS quasar positions and the FIRST catalogue 



within 2 arcsec (see e.g. Ilvezic et al.l |2002| : iLu et al.l 120071 ). resulting in a sample of 516 quasars. These 516 quasars were 



Kimball fc Ivezic 



further cross-correlated between the SDSS quasar positions and the GB6 catalogue within 1 arcmin (e.g 
20081 ) . This results in a sample of 212 quasars. The 187 quasars are thus defined as FSRQs conventionally with a spectral 
index between 1.4 and 4.85 GHz a < 0.5. After excluding two FSRQs due to the weakness of emission lines, our final sample 
consists of 185 FSRQs. The source redshift ranges from ~ 0.1 to ~ 4.0, however only one source (SDSS J081009.9-I-384756.9, 
z — 3.95) has a redshift of z > 3.0. In 79 out of 185 source, the red shift is 2: < 0.8, enabling us to measure th e broad H/3 



line, which is commonly used to estimate the black hole mass (e.g. iKaspi et al 



I2OO0I : IGu. Cao fc JianelboOll ). Based on 
the identified infrared counterpaj ts labelled in SDSS D R3 catalogue, the NIR {J,H,Ks) data are archived from the Two 



Micron All Sky Survey (2MASS, ISkrutskie et"al]l2006l ) for 98 sources among 185 sources These sources were consistently 



re-identified through cross-correlating between 2MASS and SDSS within a matched radius of 2 arcsec, which is much larger 
than the sub-arcsec positional accuracy of the SDSS and the 2MASS surveys. M oreover, we collected the Far- and near-UV 
magnitudes from the Galaxy Evolution Explorer (GALEX; Martin et al. 20051 ) Data Release 4 matched within 3 arcsec of 
SDSS positions for 117/185, and 146/185 sources, respectively. 

Our sample is listed in Table[T]and Tabled of which: (1) Source SDSS name; (2) redshift; (3) black hole mass (§ 3.2); (4) 
synchrotron peak luminosity (§ 4.4); (5) broad line region luminosity (§ 3.2); (6) bolometric luminosity (§ 3.2); (7) synchrotron 
peak frequency (§ 4.4); (8) - (10) spectral indices between the rest-frame frequencies of 5 GHz, 5000 A and 1 keV (§ 5.1). 



3 PARAMETERS DERIVATION 



3.1 Spectral analysis 



The spectra of quasars are characterized by a featureless continuum and various of broad and narrow emission lines (jVanden Berk et al 
2OOII ). For our FSRQs sample, we ignore the host galaxy contribution to the spectrum, since only very little, if any, 
starlight is observed. In a firs t step, the SDSS spectra were corrected for the Galactic extinction using the reddening map of 
Schlegel. Finkbeiner fc David ( 19981 ) and then shifted to their rest wavelength, adopting the redshift from the header of each 
SDSS spectrum. In order to reliably measure line parameters, we choose those wavelength ranges as pseudo-continua, which 
are not afi^ected by prominent emission lines, and then decompose the spectra into the following three components: 

1. A power-law continuum to describe the emission from the active nucleus. The 15 line-free spectral regions were firstly 
selected from SDSS spectra covering 1140A to 7180A for our sample, namely, 1140A - 1150A, 1275A - 1280A, 1320A - 1330A, 
I455A - 1470A, 1690A - I700A, 2I6OA - 2I8OA, 2225A - 2250A, 3OIOA - 3 C)40A, 3240A - 3270A, 3790A - 3810A, 4210A 
- 4230A, 5O8OA - 5IOOA, 5600A - 5630A, 5970A - 6OOOA, 7160A - 7180A (jVanden Berk et aLlEoOll : iForster et~aLll200ll ). 
Depending on the source redshift, the spectrum of individual quasar only covers five to eight spectral regions, from which 
the initial power-law are obtained for each source. We found that a single power-law can not give a satisfied fit for some 
low-redshift sources whose s pectra cover Ha and H/3 region. In these cases, a double power-law was adopted to obtain the 
initial power-law continuum IjVanden Berk et al.ll200ll ) with the break point around 5100A. 
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2. An Fe II template. T h e spe ctra of our sample covers UV and optical regions, therefore, we adopt the UV Fe II template 
from IVestergaard fc Wilked (|200ll '). and optical one from Veron-Cetty et al. (2004). The Fe II template obtained by Veron- 
Cetty et al. (2004) covers the wavelengt hs between 3535 and 7534 A, extending farther to both the blue and red wavelength 
ranges than the Fe II template used in iBoroson fc GreenI (|l992l ). This makes it more advantageous in modeling the Fe II 
emission in the SDSS spectra. For the sources with both UV Fe II and optical Fe II l ines prominent in the spectra, we connect 
the UV and optical templates into one template covering the whole spectra (see also Hu et al. 2008h . In the fitting, we assume 
that Fe II has the same profile as the relevant broad lines, i.e. the Fe II line width usually was fixed to the line width of broad 
II/3 or Mg II or C IV, which in most cases gave a satisfied fit. In some special cases, a free-varying line width was adopted in 
the fitting to get better fits. 

3. A Balmer con tinuum generated in the same way as Dietrich et al. (2002) (see also Hu et al. 20081 ). Grandi (1982) and 
Dietrich et al. I l|2002l ) proposed that a partially optically thick cloud with a uniform temperature could produce the Balmer 



continuum, which can be expressed as: 

^Bac ^ FbeBx{T,)[1 - e""^); (A < Abe) 

where Fbe is a normalized coefficient for the flux at the Balmer edge (Abe ~ 3646A), B\{Te 
electron temperature Te, and ta is the optical depth at A and is expressed as: 

A . 



(1) 

is the Planck function at an 



Ta — tbe{ 



Abe 



(2) 

where tbe is the optical depth at the Balmer edge. There are two free parameters, Fbe and tbe- Following Dietrich et al.l 
I 2002ll . we adopt the electron temperature to be Te = 15, 000 K. 

The modeling of above three components is performed by minimizing the in the fitting process. The final multi- 
component fit is then subtracted from the observed spectrum. The examples of the fitted power-law, Fe II lines, Balmer 
continuum and the residual spectra for the sources with low, middle and high redshift are shown in Fig. [1] The Fe II fitting 
windows ar e selected as the regi ons with prominen t Fe II line emission while no other strong emission lines, according to 
Ivestereaard fc Wilkes! (2001) and Kim et al. 1 20061 ) . The fitting window around the Balmer edge (3625 — 3645 A) is used to 
measure the contribution of Balmer continuum, which extends to Mg II line region. The Balmer continuum is not considered 
when 3625 — 3645A is out of the spectrum. Therefore, the inclusion of Balmer continuum depends on the source redshift, 
which is illustrated in Fig. [1] While the Balmer continuum should be included in the fitting for middle redshift sources, it can 
be ignored for low and high redshift sources. 

The broad emission lines were measured from the continuum subtracted spectra. We mainly focused on several prominent 
emission lines, i.e. Ly a, Ha, H/3, Mg II, C IV. Generally, two gaussian components were adopted to fit each of these lines, 
indicating the broad and narrow line components, respectively. The blended narrow lines, e.g. [O HI] AA4959, 5007A and 
[He II] A4686A blending with H/3, and [S II] AA6716, 6730A, [N II] AA6548, 6583A and [O I]A6300A blending with Ha, were 
included as one gaussian component for each line at the fixed line wavelength. The minimization method was used in fits. 
The line width FWHM, line flux of broad Ly a, Ha, H/3, Mg II and C IV lines were obtained from the final fits for our sample. 
The examples of the fitting are shown in Fig. [1] for the sources in the different redshift. 



3.2 A/bh and Lblr 

There are various empirical relations between the radius of broad line region (BLR) and the continuum luminosity, which 
can be used to calculate the black hole mass in combination with the line width FWHM of broad emission lines. However, 
there are defects when using the continuum luminosity to estimate the BLR radius for blazars since the continuum flux of 
blazars are usually doppler boosted due to the fact that the relativistic jet is oriented close to the line of sight. Alternatively, 
the broad line emission can be a good indicator of thermal emission from accretion process. Therefore for our FSRQs sample, 
we estimate the black hole mass by using the empirical relation based on the luminosity and FWHM of broad emission lines. 
According to the source redshift, we use various relations to estimate the black hole mass: Greene et al. (2005) relation for 
bro ad Ha line; Ve s terga ard et al. (2006) for broad H/3 and Kong et al. (2006) for broad Mg II and C IV lines. 



Greene fc Hcl l|2005l ) provided a formula to estimate the black hole mass using the line width and luminosity of broad Ha 



alone, which is expressed as 



Af m..)-(-n+''-^)-.m^( ^(^"^ n0.55±o.o2 . FWHM(Ha) N2.o6±o.06 
AfBH(Ha) - (2.0_o.3) X 10 ( 1342 ergs- J I 103kms-i ) ^© 



F or the sources with available FWHM and luminosity of broad H/3, the method to calculate A/bh is given bv lVestergaard fc Peterson 
1 20061) : 



MMm - 4.68 X 10« r,#^^V " f ?=^»V Mq (4) 
^ ' V 10 ergs- V VlOOOkms-i/ ^ 



In addition, Kong et al. ( 20061 ) presented the empirical formula to obtain the black hole mass using broad Mg II and 
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C IV for high redshift sources as follows, 



A/BH(MgII) = 2.9xlO« 

ylO^-^ergs-i 

Mbh(C IV) =4.6x10^ (t^^^ 
\ 10*^ erg s ^ 



0.57±0.12 /pYVHM(Mg 11)^ ^ 



V 1000 km f 



M, 



0' 



/ FWHM(C IV) 



M(7 



© 



(5) 



(6) 



1000 kms-' 

In the redshift range of our sample, Mbh can be estimated using two of above relations for 122 out of 185 FSRQs. In 
most sources (113/122), two Mbh values are consistent with each other within a factor of three. For low-redshift sources, we 
first selected the black hole mass estimated from broad H/3 line, whic h is commonly used to estimate the black hole mass 
for low-redshift sources (e.g. iKaspi et al.ll2000l : IGu. Cao fc Jiang II2OOII ). In the case that the spectral quality or the spectral 
fitting of Hz? region is poor, we instead use broad Ha line to estimate the black hole mass. Moreover, we adopted the average 
value of AfBH (Mg II) and A/bh (C IV) when both values are av ailable for one source. 



In this work, the BLR luminosity Lblr is derived following ICelotti. Padovani fc Ghisellinil (|l997l ) by scaling the strong 
broad emission lines Ly a. Ha, H/3, Mg II and C IV to the quasar template spectrum of Francis et al. (1991), in which Lya 
is used as a reference of 100. By adding the contribution of Ha with a value of 77, the total relative BLR flux is 555.77, of 
which Lya is 100, Ha 77, H/3 22, Mg II 34, and C IV 63 (|Celotti. Padovan i fc Ghisellinil ll997l : iFrancis et al.lll99lh . From the 
BLR luminosity, we estimate the bolometric luminosity as Lboi = IOLblr ()Netzeriil99Q '). 



4 THE THERMAL EMISSION 

In addition to the relativistically beamed, non-thermal jet emission, the thermal emission from the accretion disk and the 
host galaxy are expected to be present in rad io quasars. In some cases, the thermal emission can be dominated over the 
non thermal jet e mission (e.g. Landt et al. 20081 ) . We estimated the contribution of thermal emission in SEDs as follows (see 
also lLandt et al.i.2008 ). and our sample is thus refined to include only the sources with SEDs dominated by the nonthermal 
jet emission. 



4.1 The accretion disk 

Following iD'Elia. Padovani fc Landtl (|2003l ') and lLandt et al.l (|2008l ') . we calculated accretion disk spectra assuming a steady 
geometrically thin, optically thick accretion disk. In this case the emitted flux is independent of viscosity, and each element 
of the disk face radiates roughly as a blackbody with a characteristic temperature, which depends only on the mass of the 
black hole, A^bh, the accretion rate, M, and the radius of the innermost stable orbit (e.g., Peterson 1997; Frank et al. 2002). 
We have adopted the Schwarzschild geometry (nonrotating black hole), and for this the innermost stable orbit is at nn — 6rg, 
where rg is the gravitational radius defined as rg = GMbh/c^ , G is the gravitational constant, and c is the speed of light. 
Furthermore, we have assumed that the disk is viewed face-on. The accretion disk spectrum is fully constrained by the two 
quantities, accretion rate and mass of the black hole. We have calculated the accretion rate using the relations Lboi = eMc? , 
where e is the efficiency for converting matter to energy, with e ~ 6% in the case of a Schwarzschild black hole. The bolometric 
luminosity is estimated as Lboi = /~^Lblr with / the BLR covering factor, which is not well known, and can be in the 
range of ~ 5% — 30% (Maiolino et al. 2001 and references therein). As in Section [3.21 we adopt a canonical value of / ~ 10% 
(Peterson 1997). The contribution of accretion disk thermal emission are estimated by calculating the fraction of the thermal 
emission to the SED data at SDSS optical and GALEX UV region. Tentatively, we simply use a marginal value of 50% at 
most of SED wavebands to divide the FSRQs into thermal-dominated (> 50%) and nonthermal-dominated (< 50%), and 
found that the thermal emission can be dominant in 100/185, 35/185, 2/185 FSRQs for / = 5%, 10%, and 30%, respectively. 



4.2 The host galaxy 

Usu ally, the host galaxie s of radio quasars are bright ellipticals, and their luminosity only spread in a relatively narrow range 
McLure et 81112 004). We estimated the contribution from host galaxy thermal e mission using the e lliptical galaxy template 



(e. 
of 



Mannucci et al 



(|2001 ). which extends from near-IR to UV frequencies (see also Landt et al. 20081 ) . Slightly different from 



Landt et al.l lj2008l ). we use the bulge absolute luminosity in R-band estimated from A/bh — Mb. relation of iMcLure et al 
l|2004l) in combination with the estimted black hole mass in Section [3.21 

log A^bh/Mq = -0.50(±0.02)AfR - 2.74(±0.48) (7) 

The contribution of host galaxy thermal emission are estimated by comparing the calculated thermal emission with the SED 
data at SDSS optical, GALEX UV and 2MASS NIR region. The value of 50% is used to distinguish thermal-dominated with 
thermal emission > 50% and non-thermal dominated with thermal emission < 50%. We found that the thermal emission can 
be dominant in about 16 of 185 FSRQs. 
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4.3 Sample refinement 

Although the thermal emission may dominate in only small fraction of sources, i.e. accretion disk thermal emission in 35 
sources, and host galaxies emission in 16 sources, we combine the contribution of accretion disk and host galaxy to maximize 
the source number, of which the non-thermal jet emission is not dominated in SED. We simply add the expected thermal 
emission from accretion disk (using a canonical value / = 10%) and host galaxy, then compare with the SED data at SDSS 
optical, GALEX UV and 2MASS NIR region. Finally, we found the combined thermal emission can dominate over (> 50%) 
nonthermal jet emission in about 67 sources, which are then called thermal-dominated FSRQs in this work and listed in Table 
(2] The remaining 118 sources are recognized as nonthermal jet-dominated FSRQs in this work (see Table [TJ. 



4.4 



and uLv, 



peak 



The SED of each quasar was constructed from multi-band data, which covers radio (1.4 and 4.85 GHz), optical (5-8 line-free 
windows selecte d from SPSS spectra), and X-ray (IkeV) da ta. The optical continuum were picked out with five to eight 
line-free regions (jForster et al.ll200ll : IVanden Berk et al.ll2001l ) from the SDSS spectra in the source rest frame, of which no or 
only very weak Fe lines are present (see Section 3.1). The radio flux at 1.4 and 4.85 GHz are K-corrected to the source rest 
frame using the spectral index betwee n these two frequ encies. The X-ray 1 keV luminosity are obtained af ter K-correction 
and c orrecting the Galactic extinction ( Shen et al. 20061 ). The IR (J, H, and Ks) data collected from 2MASS ( Skrutskie et al. 
2006h were also added to construct SEDs, which are available for 98 FSRQs. Moreover, the Far- (for 117 sources) and near-UV 



(for 146 sources) data are also added in constructing SEDs, after correcting the Galactic extinction and K-correction using a 
spectral index of 0.5. Considering the possibility that the X-ray emission of FSRQs can be from the inverse Com pton process, 



we fit ted the data points for each source (in a v versus vh^ diagram) with a third-degree polynomial following iFossati et al. 



1 1998h . which yields an upturn allowing for X-ray data-points that do not lie on the direct extrapolation from the lower energy 
spectrum (see examples in Fig. [2]). Through fitting, we obtained the synchrotron peak frequency and the corresponding peak 
luminosity for each source. In following analysis, we will present the analysis only for the non-thermal jet-dominated FSRQs 
(see Table m. 



5 RESULTS 



5.1 



Oox plane 



The broadband properties of our sources can be firstly studied by deriving their Qox, ctro, and arx values, which are the usual 
rest-frame effective spectral indices defined between 5 GHz, 5000 A, and 1 keV. The 5000 A optical continuum fiux density is 
derived (or extrapolated when 5000 A is out of the SDSS spectral region) from the direct power-law fit on the 5-8 selected 
line-free SDSS spectral windows in the source rest frame (see Section 3.1). The 5 GHz fiux density have been k-corrected 
using the spectr al index b e tween FIRST 1.4 GHz and GB6 4.85 GHz. In Fig. [S] we present Qro - "ox relation for the sample. 
Following .Padovani et al. ( 20031 ) . three lines are indicated: Orx = 0.85, typical of 1 Jy FSRQs and LBLs; Oi-x = 0.78, the 
dividing li ne between HBLs and LBLs; and Orx = 0.70, typical of RGB BL Lac objects. Moreover, the 'HBL' and 'LBL' boxes 
defined in Padovani et al. I l|2003l) are also indicated in the figure, which repres ent the regions within 2a from the mean Oro, 
Qox, and Orx values of HBLs and LBLs in the multifrequency AGN catalog of iPadovani et al.l (|1997l ). respectively (see Fig. 
1 in Padovani et al. 2003). This "HBL box" is expected to be populated by high-energy peaked blazars, both BL Lacs and 
FSRQs. Among total 118 FSRQs, 48 sources have Qrx < 0.78, of which 28 sources locate in HBL box. In contrast, 59 sources 
among 70 Qrx > 0.78 sources are in LBL box. 



5.2 



and vLv, 



peak 



The relation of the synchrotron peak frequency and peak luminosity is presented in Fig. [l] We found only a weak anti- 
correlation between the synchrotron peak frequency i/pcak and the peak luminosity uLvpca.k with the Spearman correlation 
coefficient r = —0.161 at ~ 92% confidence level. The Upcak distribution ranges 10^^" and lO^*^-^ Hz for whole sample, and 
between 10^^ and 10^^'^ Hz for most (111/118) of sources, with (log //peak) = 14.41 ± 0.74 Hz for whole samp le. We found 



I'peak > 10^^ Hz in 25 sources (see Fig. [Jjl, which is larger than the typical value of FSRQs I Fossati et al. 19981'). As outliers 



to bl azar sequence, the blue quasars are supposed to have large /^pcak and high i/Li/poak as well (see e.g. iPadovani et al 
2003[ ). To further understand the nature of these high !^pcak FSRQs, we combine our sample with the sample of 170 BL Lac 



objects in Wu. Gu fc Jiang 1 20081 ). A significant anti-correlation between z/pcak and i/Lfpeak is present with the Spearman 
correlation coefficient r = —0.343 at ^ 99.99% confidence level for the combined sample of 288 blazars, which is shown 
in Fig. [5] This anti-correlation covers about seven order of magnitude in i^peak for the combined sample. Our FSRQs with 
high Upcak > W^'' Hz have systematically higher peak luminosity than that of BL Lacs, with (log v'Lvpcaia) ~ 46.41 ± 0.94 
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compared to 44.90 ± 0.71 for BL Lacs at same Vpcak range, and the peak luminosity of some FSRQs are out of the range 
covered by BL Lacs. However, we found that the thermal emission from accretion disk can be dominated over nonthermal 
jet emission in 11 of 25 sources with t/peak > lO^'^ Hz if a BLR covering factor of 5% is assumed. Therefore, the possibility 
that the contribution from the thermal emission causes the high synchrotron peak frequency can not be completely excluded. 
Although high t'peak are found in our FSRQs, they do not reach the extreme v alue of HBLs, which is also found in DXRBS 
sample when comparing high i/pcak FSRQs with BL Lacs ( Padovani et al. 20031 ). At the lower-left corner of Fig. [S] there are 
some FSRQs with relatively low synchrotron peak frequency as well as low lumi nosity, comparable to so me low luminosity 
LBLs. While only a weak anti-correlation present between I'pcak and uLvpak for IWu. Gu fc Jiand l|2008f ) BL Lacs sample, 
it becomes significant when combinin g with our sampl e. To some extent, this reflects the importance of sample selection in 
investigating the blazar sequence (e.g. iPadovani Il2007^ . Although the anti-correlation is significant in Fig. [S] we notice that 
the scatter is significant, mainly caused by the sources in the lower-left corner. However, the extreme sources with high z^pcak 
as well as high peak luminosity (i.e. at upper-right corner) are still lacking. On the other hand, the large scatter may imply 
that other parameters may be at work, and the scatter can be much lower once this parameter is included. 

We investigate the relationship between arx and Vpcak for our sample in Fig. [G] A significant anti-correlation is found with 
a Spearman correlation coefficient r = —0.443 at confidence level S> 99.99%. However, the considerable scatter is also present 
with the scatter in z^peak more than one order of magnitude for any given value of Ofrx- The mean I'peak of FSRQs in HBL box 
is only slightly larger than that of sources out of the HBL box with a factor of four. Consequently, the FSRQs in HBL box do 
not exclusively have high t'peak. The i^peak of souces in HBL box are comparable to those with < 0.78 but out of HBL box. 
This actually can be refiected from the relationship between a^o and I'peak shown in Fig. [T] The significant anti-correlation is 
found at a confidence level of ~ 99.9%, however the scatter is significant. At any given value of ctro, the source could have a 
low Vpcak or a high i/pcak- As a result, the sources out of HBL box with Qrx < 0.78 are not necessary to have lower or higher 
Vpeak than that of sources in HBL box, though the latter show a wider Qro range. Consistent with the anti-correlation in Fig. 
[6l the mean Vpcak value of FSRQs with Orx < 0.78 is larger with a factor of four than that of Orx > 0.78 sources. Owing to 
the large scatters in both Qrx — i-'peak and Qro — Vpcak relations, it seems that neither solely Qrx nor Qrx - Qrx combination can 
precisely predict the Vpcak value for our present FSRQs sample. 



5.3 (Upcak, uLUpcak) & (A^bh, ibol/i/Edd) 

The estimation of black hole mass and BLR luminosity (then the corresponding bolometric luminosity) enable us to investigate 
the relationship between the synchrotron emission and the accretion process. In Figs. [8] and |9] we show the relation of Vpcak & 
Mbh, and Vpcuk & iboi/^Edd, respectively. The black hole mass ranges from 10^ * to 10^"'^ Mq with most of sources (~ 86%; 
102 out of 118 sources) in the range of 10*'^ — 10^" Mq, while the eddington ratio I/boi/^'Edd ranges from ~ 10~^ to ~ 10°'®, 
with most of source (~ 94%; 111 out of 118 sources) in the range of 0.01 to 1. No strong correlations are found either between 
and Mbh, or between z/poak & iboi/^Edd. 
We found a significant correlation between the black hole mass Mbh and the synchrotron peak luminosity z/Li^pcak with 
the Spearman correlation coefficient r = 0.724 at ~ 99.99% confidence level, which is shown in Fig. [lOl Moreover, a significant 
correlation between lyLi/pcuk and Lboi/iEdd is also found with the Spearman correlation coefficient r = 0.842 at S> 99.99% 
confidence level (see Fig. Ill[). Both correlations still present when we perform the Spearman partial correlation analysis to 
exclude the common dependence on redshift. The ordinary least-square (OLS) bisector linear fit to fLvpcak and Lboi/iEdd 
gives, 

log /.il/peak = (1.73 ± 0.11) log ibol/iEdd + (47.92 ± 0.11) (8) 

This result indicates that the jet physics may be tightly related with the accretion process. 

We plot the ratio of BLR luminosity to synchrotron peak luminosity I/BLR/fii^poak and synchrotron peak frequency 
t'pcak in Fig. 1121 from which a signiffcant correlation between these two parameters is found with the Spearman correlation 
coefficient r = 0.480 at ^ 99.99% confidence level. However, there is no strong correlation between the synchrotron peak 
frequency and BLR luminosity. 

The jet-disk relation can be further explored through the relationship between i/Li/pcaii and Lblr shown in Fig. 1131 A 
signiffcant correlation is present with a Spearman correlation coefffcient of r = 0.909 at confidence level ^ 99.99%, which 
remains in partial correlation analysis to exclude the common dependence on redshift. The ordinary least-square (OLS) 
bisector linear fit gives, 

log z/ii/peak = (0.95 ± 0.06) log Lblr + (3.14 ± 2.84) (9) 

Our result implies a tight relation between jet and disk, which is consistent with previous findings in various occasions, from 
the strong correlations either between the radio emission and emission line luminosity (e.g. Rawlings et al. 1989; Cao & Jiang 
2001), or between the emission line luminosity and jet kinetic power in different scales (Rawlings & Saunders 1991; Celotti 
& Fabian 1993; Wang et al. 2004; Gu, Cao & Jiang 2009). Moreover, our result is consistent with the nearly proportional 
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relation, Q oc L^'^^''^ , found between the jet bulk kinetic power and narrow line luminosity (Rawlings & Saunders 1991). 
However, it should be noted that the synchrotron peak luminosity of blazars are usually Doppler boosted, therefore, it may 
not be a good indicator of jet power. It would be necessary to improve our result using intrinsic parameters, either the intrinsic 
synchrotron peak luminosity after eliminating the beaming effect, or the jet power. 



6 DISCUSSION 



According to the blazar sequence, FSRQs with high synchrotron peak frequency, e.g. t-pcak > 10^® Hz, and X-rays dominated 
by synchrotron emission are not expected to exist. Due to its importance, such objects have been extensively searched (see 
PadovatiilboO? , and references therein) . The discoveries of such blaza rs have been claimed, however, mainly on the basis of their 
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broad spectral indices, i.e. the ratios of radio/optical/X-ray fluxes IjPadovani et al.ll2002l : iBassani et al. 
20071 ). In this sense, the X-ray spectroscopy is required to confirm the nature of the X-ray emission in these blazars. In 
its simple st form, a X-ra y spec tral index of a < 1 indicates a inverse Compton origin, while a > 1 for synchrotron X-ray 



Maraschi et al. 1 2008! ) investigate the X-ray spectra for a sample of 10 X-ray selected FSRQs from the Einstein 



Medium Sensitivity Survey (EMSS) and four controversial sources claimed to have synchrotron X-ray emission. They found 
that, in the case of the EMSS broad line blazars, X-ray selection does not lead to find sources with synchrotron peaks in the 
UV/X-ray range, as was the case for X-ray-selected BL Lacs. Instead, for a wide range of radio powers all the sources with 
broad emission lines show similar SEDs, with synchrotron components peaking below the optical/UV range. Moreover, the 
authors argued that four 'anomalous' blazars are no longer 'anomalous' after a complete analysis of Swift and INTEGRAL 
data, with two sources having inverse Compton X-rays, one source being HBL, and the remaining one being narrow line 
Seyfert 1 galaxy without unambiguous evidence of X-ray emissio n from a rela tivistic jet. Similarly, the XMM — Newton and 
Chandra X-ray spectroscopy of 10 FSRQs were investigated bv lLandt et al. 1 lioos), which are candidates to have an X-ray 
spectrum dominated by jet synchrotron emission. However, the authors failed to find FSRQs with X-ray spectra dominated 
by jet synchrotron emission, instead, the X-rays are either from inverse Compton or are at transition between the synchrotron 
and inverse Compton jet components as in IBLs. So, despite the efforts to search for objects which may violate the sequence 
trends no strong outliers have been found (Maraschi et al. 2008). As a severe challenge to blazar sequence, the selection of 
important candidates of high !/peak luminous FSRQs ar e refined to choose highly core-dominated radio quasars with low radio 
core to X-ray luminosity ratios, e.g. log (icoro/ix) 5 ( Landt et al. 20081 ) . The inverse Compton emission is expected to peak 
at 7-ray frequenci es, therefore, the high energy-peaked FSRQs could be prime targe ts for the Fermi Gamma-ray Telescope 
I Landt et alJ 20o"3 l . Although the high energy-peaked FSRQs are not firmly found yet, Ghisellini fc Tavecchiol (2003) proposed 
a model to explain the existence of blue quasars. In their scenario, the jet dissipation region is out of the broad line region, 
resulting in a much reduced energy density of BLR photons in jet region. Therefore, the cooling due to the inverse Compton 
process is not severe, causing a high synchrotron peak frequency though the source luminosity is high. 

Although the anti-correlation between Vpcsk and //Li^peak for the combined sample of 288 blazars (Fig. [5} is significant, it 
largely differs with the blazar sequence in its significant scatter, which indicates that the synchrotron peak luminosity can not 
explicitly determine the synchrotron peak frequency, and vice verse. The advantages of using the synchrotron peak luminosity 
is that the most of synchrotron emission are radiated at synchrotron peak frequency, at which the luminosity can be a good 
indicator of synchrotron emission. Since the synchrotron peak frequency varies from source to source, the luminosity at fixed 
wavebands (e.g. optical) is actually from the different portion of source SED. The defect of the synchrotron peak luminosity 
lies in the contamination from beaming effect, which precludes it to well indicate the intrinsic source power. Only when the 
Doppler boosting is known for each source, the intrinsic source power can be obtained. This can not be performed at present 
stage. However, this effect has been explored through eliminating the Doppler boosting in blazar samples. Interestingly, 
recent studies have shown that the negative correlation between Vpcsk and i^Lt'pcak is likely an artefact of Doppler boosting 
I Nieppola et al. 20081 : Wu. Gu fc Jians 20081 ). According to these authors, the negative correlation is not present when the 
intrinsic parameters are used, conversely, a positive correlation is claimed. The key point in these studies is the strong anti- 
correlation between the Doppler factor and synchrotron peak frequency in the way that the sources with their synchrotron 
peak at low energies are sig nificantly more boosted than high //peak sources, whic h is found either from t he variability Doppler 
factor ( Nieppola et al. 2008 ) or from the one estimated with empirical relation ( Wu. Gu fc Jiang 20081 ) . 

The synchrotron peak frequency Up^sA. oc BS^'^^^^^, where B is the magnetic field, 5 the Doppler factor, and 7peak a 
characteristic electron energy that is determined by a competition between accelerating and cooling processes. Different from 
BL Lacs, the external inverse Compton scattering is thought to be the dominant cooling process in FSRQs, especially that 
upon BLR photons. Through model fitting to blazar SEDs, the electron peak energy is found to be well anti-correlated with 
the total energy density (radiative and magnetic), which is thought t o be the phy sics behind the phenomenological blazar 
sequence ( Ghisellini et al. 19981 : Ghisellini. Celotti fc Costamante 2002 ; Celotti fc Ghisellini 2008f l . More energy density inside 
the jet cause a more severe cooling, resulting in a smaller 7peak then I'peak- According to the equations we used to estimate the 
black hole mass, we have the BLR radius approximately with -Rblr oc I/blr for Hq, H/3, Mgll, CIV broad lines. Consequently, 
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the energy density of BLR photons iiblr is expected to be proportional to L^^'^ ■ We expect to see an anti-correlation between 
J'pcak and the BLR luminosity Lblr- However, we failed to find any correlation between Vpcak and Lblr- Several factors may 
erase the expected anti-correlation, e.g. the scatters in the derivation of BLR luminosity from individual lines, the accuracy of 
empirical relation in estimating BLR radius, and the inclusion of B and S (vary from source to source) in Upca.k- As a result, 
the positive correlation between z^pcak and LBLR/i^Li/pcak can be partly (if not all) the result of the weak anti-correlation 
between Vpe^k and vLupeaii (see Fig. I12|l . Nevertheless, it indicates that FSRQs with higher ratio of disk to jet emission could 
have higher peak frequency. 

The FSRQs in our sample are defined from the flat spectrum between 1.4 and 4.85 GHz a < 0.5. Howeve r, this definitio n 
may influenced by several factors. It is well known that FSRQs usually show strong and rapid variability (e.g. Gu et al. 20061 ) . 
Another factor is the different resolution at 1.4 and 4.85 GHz. FIRST 1.4 GHz data are obtained from VLA observations, which 
have much higher resolution than Green Bank telescope observations for 4.85 GHz data. In these respects, the simultaneous 
multi-band observations with same telescope conflguration (same resolution) is required to calculate the radio spectral index, 
and then to understand the nature of sources. 



7 SUMMARY 

We have constructed a sample of 185 Flat Spectrum Radio Quasars (FSRQs) by cross-correlating the lShen et al. I l|2006l ) SDSS 
DR3 X-ray quasar sample with FIRST and GB6 radio catalogues. From the spectrum energy distraction (SED) constructed 
using multi-band (radio, optical. Infrared and X-ray) data, we derived the synchrotron peak frequency and peak luminosity. 
The black hole mass Mbh and the broad line region (BLR) luminosity (then the bolometric luminosity Lboi) were obtained 
by measuring the line-width and strength of broad emission lines from SDSS spectra. We define a subsample of 118 FSRQs, 
of which the nonthermal jet emission are thought to be dominated over thermal ones from accretion disk and host galaxy. 
The various correlations were explored for this subsample. The main results are summarized below. 

1. A weak anti-corre lation is found be t ween the synchrotron peak frequency and peak luminosity. When combining our 
FSRQs sample with the Wu. Gu fc Jian j 1 2008h sample of 170 BL Lac objects, a significant anti-correlation between the 
synchrotron peak frequency and luminosity apparently presents covering about seven order of magnitude in t'poak. However, 
the anti-correlation differs with the blazar sequence in the large scatter. 

2. We found 25 FSRQs having synchrotron peak frequency Vpeak > 10^^ Hz, which is higher than the typical value for 
FSRQs. These sources with high z/poak could be the targets for the Fermi Gamma-ray telescope. At similar peak frequency, 
the peak luminosity of FSRQs with i^peak > 10^^ Hz is systematically higher than that of BL Lac objects, with some FSRQs 
out of the range covered by BL Lac objects. Though high z^pcak are found in some FSRQs, they do not reach the extreme 
value of BL Lacs. 

3. No strong correlations are found either between the synchrotron peak frequency and black hole mass, or between the 
synchrotron peak frequency and the Eddington ratio. The peak luminosity is found to be tightly correlated with both black 
hole mass and the Eddington ratio indicating that the jet physics may be tightly related with the accretion process, which is 
further confirmed by the tight correlation between the synchrotron peak luminosity and the BLR luminosity. 
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Table 1. Nonthermal jet-dominated FSRQs. (1) Source SDSS name; (2) redshift; (3) black hole mass; (4) synchrotron peak luminosity; 
(5)broad line region luminosity; (6) bolometric luminosity; (7) synchrotron peak frequency; spectral indices between the rest-frame 
frequencies of (8) 5 GHz and 5000 A, (9) 5000 A and 1 keV, (10) 5 GHz and 1 keV. 
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46.50 


14.11 
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1.320 


0.815 


J130603.4-f552943.7 


1.601 


9.53 


46.90 


46.25 


47.25 


14.38 


0.473 


1.489 


0.818 
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SDSS name 


z 


Log Mbh 
(■^0) 


Log J^L^poak 

(ergs ^) 


Log Lblr 
(ergs ^) 


Log Lboi 
(ergs ^) 


Log J/peak 

(Hz) 




Oox 


Qrx 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


v) 


(8) 


(9) 


(10) 


J130909.8+555738.5 


1.629 


9.64 
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46.31 
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14.80 


0.450 


1.665 
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0.669 
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44.59 


45.59 
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0.879 


J133437.4+563147.9 
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14.33 


0.565 


1.640 
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9.56 
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0.933 
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46.13 


14.70 
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0.829 
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9.11 
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14.86 
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1.271 


0.813 


J135817.6+575204.4 


1.373 


9.72 
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13.92 


0.210 


1.693 


0.713 


J140126. 2+520834. 6 


2.972 


9.69 


47.86 


47.19 


48.19 


15.01 


0.219 


1.626 


0.696 


J141159.7+423950.4 


0.886 


9.60 


46.25 


45.72 
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14.85 


0.382 


1.629 


0.805 


J141324. 2+530525. 7 


0.456 


8.10 


44.57 


43.55 
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1.289 
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1.559 


9.68 
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0.914 


J142020.7+462441.0 


1.245 


9.50 
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45.30 
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1.158 


9.48 


46.08 
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0.392 
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16.27 
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0.591 
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1.701 


9.61 
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1.101 


9.16 
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13.93 
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14.36 
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0.771 
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9.68 
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9.28 
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14.71 


0.368 


1.301 
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1.941 


9.92 
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13.12 


0.610 


1.728 


0.989 
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Table 2. Thermal-dominated FSRQs. (1) Source SDSS name; (2) redshift; (3) black hole mass; (4) synchrotron peak luminosity; (5)broad 
line region luminosity; (6) bolometric luminosity; (7) synchrotron peak frequency; spectral indices between the rest-frame frequencies of 
(8) 5 GHz and 5000 A, (9) 5000 A and 1 keV, (10) 5 GHz and 1 keV. 
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J152556.2+591659.2 


0.955 


9.84 


46.19 


45.66 


46.66 


13.91 


0.461 


1.485 


0.808 


J152942.2+350851.5 


0.287 


8.73 


44.62 


44.06 


45.06 


14.01 


0.506 


1.283 


0.769 


J153404.4+580059.4 


1.572 


9.88 


46.73 


46.44 


47.44 


13.83 


0.274 


1.538 


0.702 


J154232.0+493842.7 


0.590 


8.96 


44.87 


44.65 


45.65 


14.47 


0.534 


1.183 


0.754 


J155853.6+425817.6 


0.874 


9.60 


45.89 


45.38 


46.38 


13.88 


0.452 


1.396 


0.772 


J161035.4-I-480022.3 


0.247 


8.26 


44.33 


43.42 


44.42 


13.68 


0.413 


1.313 


0.718 


J161655.6-f362134.4 


2.265 


9.90 


47.23 


46.91 


47.91 


14.64 


0.473 


1.301 


0.754 


J163856.5-f433513.0 


0.339 


8.97 


44.81 


44.04 


45.04 


15.09 


0.472 


1.417 


0.792 
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Table 2. Continued. . . 



SDSS name 


z 


Log Mbh 


Log J^L^poak 


Log Lblr 


Log Lboi 


Log J/peak 






Qrx 






(Mo) 


(ergs-^) 


(ergs-^) 


(ergs~i) 


(Hz) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


J164051. 2+314329.8 


0.958 


9.50 


45.66 


45.38 


46.38 


14.16 


0.559 


1.316 


0.816 


J164147.6+393503.4 


0.539 


9.04 


45.18 


44.68 


45.68 


15.95 


0.499 


0.986 


0.664 


J164452.6+373009.3 


0.758 


9.99 


46.03 


45.32 


46.32 


13.59 


0.473 


1.533 


0.832 


J164544.7+375526.3 


0.598 


9.34 


45.45 


44.85 


45.85 


14.24 


0.464 


1.482 


0.809 


J165801.4+344328.5 


1.939 


9.85 


47.41 


46.69 


47.69 


13.71 


0.523 


1.624 


0.896 


J170112.4+353353.3 


0.501 


8.70 


44.62 


44.42 


45.42 


14.95 


0.574 


0.722 


0.624 


J172051. 2+620944.6 


1.010 


9.65 


46.19 


45.45 


46.45 


13.46 


0.374 


1.189 


0.650 
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Figure 1. The examples of spectral analysis: (A) low redshift FSRQ SDSS J080131. 9+473616.1 (z = 0.157); (B) middle redshift FSRQ 
SDSS J074559.3+331334.5 {z = 0.610); (C) high redshift FSRQ SDSS J081303.8+254211.1 {z = 2.024). In each panel, the top and 
middle black lines are the original and the continuum-subtracted spectrum, respectively. The top dotted red line represents the power- 
law continuum. The middle blue lines are the individual line components in multi-line spectral fitting, and the middle solid red line 
is integrated lines fitting. The green spectral region in the continuum-subtracted spectrum is used to fit Fe II emission and power-law 
continuum, and the bottom green line is Fe II emission, which is shifted downwards with arbitrary unit for the sake of presentation. The 
Balmer continuum is indicated as solid blue lines in the bottom of panel (B), which began at the Balmer edge 3646A, however, it is not 
included in the case of (A) and (C) due to the wavelength coverage (see text for details). 
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Figure 2. Representative three-degree polynomial fits to the spectral energy distributions of FSRQs in our sample. Open rectangles are 
data from FIRST, GB6, SDSS and RASS. Stars represent data from 2MASS, and triangles for GALEX UV data. The solid line is the 
SED fit, while the dotted one (HG) is the calculated emission from host galaxy. The dashed line (AD) is the expected emission from the 
accretion disk at assuming / = 10% (see text for details). The synchrotron peak frequencies for the two sources are fc'peak = 1-86 x 10^^ 
Hz (SDSS J150324.7+475830.4, z = 0.345; A) and i-peak = 2.88 x lO^^ Hz (SDSS J083148.9+042939.2, z = 0.174; B). See text for details. 
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Figure 3. The Oro & ctox relation. The effective spectral indices are defined in the usual way and calculated between the rest-frame 

frequencies of 5 GHz, 5000 A, and 1 kcV. The solid circles are 118 nonthcrmal-dominatcd FSRQs, while the open ones for 67 thermal- 
dominated sources (sec text for details). The dotted lines represent, from top to bottom, the loci of a^x = 0.85, typical of 1 Jy FSRQs 
and LBLs; Qrx = 0.78, the dividing line between HBLs and LBLs; and arx = 0.70, typical of RGB BL Lac objects. The 'HBL' and 'LBL' 
boxes defined in Padovani et al. (2003) are indicated by the solid lines and marked accordingly, which represent the regions within 2a 
from the mean aro, ctox, and arx values of HBLs and LBLs, respectively (see also Fig. 1 in Padovani et al. 2003). 
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Figure 4. The synchrotron peak frequency and the peak luminosity for our sample. The symbols are same as in Fig. |3] 
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Figure 5. T he synchrotron p e ak fre quency and the peak luminosity. The circles are same as in Fig.fS] while the stars are 170 BL Lac 
objects from lWu. Gu fc Jiand ||2008|1 . 
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Figure 6. The synchrotron peak frequency and a,x- The circles represent FSRQs with Oix > 0.78, the triangles for FSRQs in HBL box, 
and the rectangles are FSRQs out of HBL box however with < 0.78. The solid symbols arc for 118 nonthermal-dominated FSRQs, 
while the open ones for 67 thermal-dominated sources. 
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Figure 7. The synchrotron peak frequency and a^o- The symbols are same as those in Fig. |6] 
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Figure 8. The synchrotron peak frequency versus the black hole mass. The symbols are same as in Fig. |3] 
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Figure 9. The synchrotron peak frequency versus ^bol/^Edd- The symbols are same as in Fig. [3] 
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Figure 10. The synchrotron peak luminosity versus the black hole mass. The symbols are same as in Fig. [3] 
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Figure 11. The synchrotron peak luminosity versus i^bol/^Edd- The symbols are same as in Fig. [3] The solid line is the OLS bisector 
linear fit for 118 nonthermal-dominated FSRQs. 
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Figure 12. The synchrotron peak frequency versus iBLn/'^ii'pcak- The symbols are same as in Fig. [3] 
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Figure 13. The vLi/peak versus Lblr- The sohd line is the OLS bisector linear fit for 118 nonthermal-dominated FSRQs. The symbols 
are same as in Fig. |3] 



